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STEVE ZELDITCH 

Abstract. We show that the expected limit distribution of the real zero set of a Gaussian 
random linear combination of eigenfunctions with frequencies from a short interval ('asymp- 
totically fixed frequency') is uniform with respect to the volume form of a compact Riemann- 
ian manifold (M, g) . We further show that the complex zero set of the analytic continuations 
of such Riemannian random waves to a Grauert tube in the complexification of M tends to 
a limit current. 
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This article is concerned with the real and complex zero sets of Riemannian random waves 
on a real analytic Riemannian manifold (M, g). To define Riemannian random waves, we fix 
an orthonormal basis of real- valued eigenfunctions of the Laplacian A g of (M,g), 

4^! A;' r -A ; . (<Px i ,<P\ k ) = 6jk, 

+-> . 

and define Gaussian ensembles of random functions / = . Cj(p\ j of the following two types: 

• The asymptotically fixed frequency ensemble TCi x , where I\ = [A, A + 1] and where 
7ii x is the vector space of linear combinations 

CO ! j:A,e[A,A+l] 

^ ! of eigenfunctions with Xj (the frequency) in an interval [A, A + 1] of fixed width. 

^ \ (Note that it is the square root of the eigenvalue of A, not the eigenvalue, which is 

O ' asymptotically fixed). 
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• The cut-off ensembles / H[o,x) where the frequency is cut-off at A: 

j:Aj<A 

By random, we mean that the coefficients Cj are independent Gaussian random variables 
with mean zero and with the variance defined so that the expected L 2 norm of / equals one. 
Equivalently, the real vector spaces W[o,a], resp. Tii x are endowed with the inner product 
(u, v) = J M uvdV g (where dV g is the volume form of (M,g)) and random means that we 
equip the vector spaces with the induced Gaussian measure. Our main results given the 
asympototic distribution of real and complex zeros of such Riemannian random waves in the 
high frequency limit A — > oo. 

The real zeros are straightforward to define. For each fx G T~i[o,\] or 7i Ix we associated to 
the zero set Zf x = {x G M : f\(x) = 0} the positive measure 

(\ZfM)= I V^- 1 , (3) 

Jz fx 
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where dTL m ~ 1 is the induced (Hausdorff) hypersurface measure. Our first result (Theorem 
[1]) shows that the normalized expected limit distribution iE|.Zyj of zeros of random Rie- 
mannian waves tends to the volume form dV g as A — > oo. The result is the same in the 
cutoff and fixed frequency ensembles, but the proof is simpler in the former. Further, in the 
fixed frequency case, if one chooses a random sequence f\ N , i.e. the elements are chosen 
independently from the frequency intervals [N, N + 1] and at random from Hwjv+i], then 

almost surely, £ J2k=i J^\ Z h k \ ~* dv r 

The complex zeros are defined by analytic continuation. Since (M, g) is assumed to be 
real analytic, it admits a complexification Mc and the eigenfunctions can be analytically 
continued to Mc- Their real zero hypersurfaces extend to complex nodal hypersurfaces 
in Mq. Our second result (Theorem [3]) determines the limit distribution of these random 
complex nodal hypersurfaces, and shows that the limit distribution is the same as the one 
determined in [Z3J for complex zeros of analytic continuations of ergodic eigenfunctions, 
e.g. eigenfunctions of A g when the geodesic flow of (M, g) is ergodic. This corroborates 
the random wave hypothesis as applied to complex nodal sets, i.e. the conjecture that 
eigenfunctions of chaotic systems resemble random waves. Note that in the real domain, 
the distribution of nodal lines of eigenfunctions in the ergodic case is very much an open 
problem. The motivation to study complex zeros of analytic continuations of eigenfunctions 
comes from the fact one has more control over complex zeros than real zeros of (deterministic) 
eigenfunctions. 

0.1. Statement of results on real zeros. We now state our results more precisely. Let us 
first consider the real zero sets of eigenfunctions of A for the standard sphere (S m ,go). Let 
TLn C L 2 (S m ) denote the real d n- dimensional inner product space of spherical harmonics of 
degree N. The eigenvalue is given by 

(X s N m r = N(N + m-l) = (N + ^r-(^ 

where (3 = is the common Morse index of the 2n periodic geodesies and 

fm + N - 1\ fm + N - 3 
dN ~{ N )-{ N-2 

We choose an orthonormal basis {(PNj}j=i for Hn- For instance, on S 2 one can choose the 
real and imaginary parts of the standard Y^'s. In the fixed frequency ensemble, we endow 
the real vector space Hn with the Gaussian probability measure 7^ defined by 

(1 \ d N /2 d x 
^) e~ d ^ 2 dc, f = Y, c 3 ( PW> d N = dimH N . (4) 

Here, dc is G?Ar-dimensional real Lebesgue measure. The normalization is chosen so that 
E 7JV (/, /) = 1, where E 7JV is the expected value with respect to 7^. Equivalently, the d^ 
real variables Cj (j = l,...,d^) are independent identically distributed (i.i.d.) random 
variables with mean and variance ; i.e., 
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We note that the Gaussian ensemble is equivalent to picking fa e Tin at random from the 
unit sphere in T~Cn with respect to the L 2 inner product. The latter description is more 
intuitive but it is technically more convenient to work with Gaussian measures. In the cutoff 
ensemble, we put the product Gaussian measure U.^ =1 j n on @^ =1 7i n . 

We now consider the analogous constructions on a general compact Riemannian manifold 
(M, g) of dimension m. As mentioned above, and as defined more precisely in $H the analogue 
of the space Hn of spherical harmonics of degree N is played by the space TCi N of linear 
combinations of eigenfunctions ([1]) with frequencies in an interval fa := [N, N + 1]. The 
precise decomposition of R into intervals is not canonical on a generic Riemannian manifold 
and the results do not depend on the choice. We choose fa = [N, N + 1] only for notational 
simplicity. In the special case of Zoll manifolds (all of whose geodesies are closed), there is a 
canonical choice (an eigenvalue cluster decomposition) which is described in §[[]). Henceforth 
we abbreviate Ti^ = Ti,j N on general Riemannian manifolds. We continue to denote by 
{fNj}j=i an orthonormal basis of Hn where = dimWjv- We equip it with the Gaussian 
measure (jl]) and again denote the expected value with respect to (HniIn) by E 7JV . 

Our main result on real Riemannian random waves is to determine the expected value 
E jN \Zf N \. It is a positive measure satisfying 

(E 7N \Z fN \,i>) = E 7N X$, (5) 

where is a 'linear statistic', i.e. the random variable 

X»(fa) = (iP,\Z fN \), ^eC(M) (6) 

considered in 

Theorem 1. Let (M,g) be a compact Riemannian manifold, let Ti.\o,x] be the cutoff ensemble 
and let (7"^v,7at) be the ensemble of Riemannian waves of asymptotically fixed frequency. 
Then in either ensemble: 

(1) For any (M,g), lim^ ±E JN {\Z fN \, ^) = f M ipdV g . 

(2) For a real analytic (M,g), Var{±X$)) < C. 

We restrict to real analytic metrics in (2) for the sake of brevity. In that case, the 
variance estimate follows easily from a result of Donnelly- Fefferman on volumes of real nodal 
hypesurfaces of real analytic (M,g). 

In the case of the standard metric go on S m , it is obvious apriori that E lN \Zf N \ = CNdV go , 
the constant Cn being the expected volume of the zero sets. The expected volume was first 
determined by P. Berard by a different method. The theorem above shows that asymptoti- 
cally the same result holds on any compact Riemannian manifold. 

A much better variance estimate for S m was obtained by J. Neuheisel in his (unpublished) 
Hopkins PhD thesis |Nej , which shows that the variance tends to zero at least at a rate 
for a certain 5 > 0. It is very likely that one could prove the same (or a better) variance 
estimate on a general C°° (M,g), but that would require a study of the pair correlation 
function of zeros which would take us too far afield from our main purpose. We plan to carry 
them out on a different occasion. Relatively sharp variance estimates for eigenfunctions on 
arithmetic tori are given in |RWt IQRWj . 

An immediate consequence, by the Kolmogorov strong law of large numbers, is a limit law 
for random sequences of random real Riemannian waves. By a random sequence, we mean 
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an element of the product probability space 

Hoo = U^ =1 Hn, 7oo = n^ =1 7Ar. (7) 

Corollary 2. Let (M,g) be a compact real analytic Riemannian manifold, and let {f N } 
be a random sequence in |?|). Then 

1 N 1 

— Y^t- \Zf n \ -> dV g almost surely w.r.t. (Hoo,^). 

It is natural to conjecture that -£-\Zf N \ — > dV g almost surely without averaging in N, but 
the proof would again require a stronger variance estimate than we currently possess. 

0.2. Statement of results on complex zeros. We now turn to results on complex zeros 
of analytic continuations of eigenfunctions. By a theorem of Bruhat- Whitney [BW], an 
analytic manifold M admits a complexification Mc into which M embeds as a totally real 
submanifold. Associated to g is a plurisubharmonic exhaustion function p(() which measures 
the square of the distance to the real subset M. The sublevel set M T = {( 6 M c : v / p(C) < T } 
is known as the Grauert tube of radius r (cf. |Grj. IGSll IGS21 ILSlj ). 

It was observed by Boutet de Monvel IBou that eigenfunctions can be analytically con- 
tinued to the maximal Grauert tube as holomorphic functions (p%. . Thus, we can complexify 
the Gaussian random waves as 

d N 

/iV = z2 C Nj<fNj- 
J'=l 

We note that the coefficients Cjvj are real and that the Gaussian measure on the coefficients 
remains the real Gaussian measure 7at. 

Our next result determines the expected limit current of complex zeros of f^. The current 
of integration over the complex zero set 

Zjc ={(eM c :f% = 0} 

is the (1,1) current defined by 



I N 



for smooth test forms of bi-degree (m — l,m — 1). In terms of scalar functions ip we may 
define Z f c as the measure, 

m— 1 



([Z f cU)= / ^V(m-l)! 



where uj 9 = iddp is the Kahler metric adapted to g. 



Theorem 3. Let (M, g) be a real analytic compact Riemannian manifold. Then for either 
of the ensembles of Theorem^ we have 

E 7iV (jj[Zfc]j - l -8d\i\ gi weakly in V'^(B* £ M). 
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As mentioned above, this result shows that the complex zeros of the random waves have 
the same expected limit distribution found in |Z3j for real analytic compact Riemannian 
manifolds with ergodic geodesic flow. 

0. 3. The key objects in the proof. The principal objects (for the asymptotically fixed 
frequency ensembles) are the two point functions 

Ui N (x,y) = E 7N (f N (x)f N (y)) = ^ <f Xj (x)<f Xj (y), (8) 

1. e. the spectral projections kernel for vA, and their analytic extensions to the totally real 
anti-diagonal in Mc x M<c defined by 

imc,c)= E iVi(or- (9) 

j:\jeh 

A key point is that the latter kernels are very much off the diagonal for non-real (, so that 
the kernels grow at an exponential rate. In the cutoff ensemble, the spectral projections 
kernels are replaced by 

K[o,N\(x,v) = ^ 1N (fN(x)f N (y)) = E VAiO^Aid/). (10) 

j:A j6 [0,iV] 

and similarly for the complexification. 

In the real domain, the distribution of zeros of random Riemannian waves is obtained 
by using formalism of [BSZ11 IBSZ2] to express the density of zeros in terms of the kernels 
Hi N (x,y). We then use the spectral asymptotics of these kernels and their derivatives to 
derive the limit distribution of zeros. A more in-depth analysis of their off-diagonal decay 
could give bounds on the variance, but as mentioned above we postpone that to a later 
occasion. In the complex domain, the spectral asymptotics have not been studied before. 
The asymptotics are more difficult than in the real domain and have an independent interest. 
On the other hand, the link between these kernels and the zero distribution is simpler, and 
we use the Poincare-Lelong method of [BSZ3J rather than [BSZ1 , BSZ2J. 

0.4. Discussion. This is the first article discussing zeros of Riemannian random waves of 
asymptotically fixed energy. We digress to compare our definitions and results on Riemann- 
ian waves to other definitions and results in the literature. 

The subject of random polynomials and Fourier series and their zeros is classical; see [PWJ 
for one of the classics. For a contemporary treatment of Gaussian random functions in a 
geometric setting, see [ATJ. Aside from its pure mathematical interest, Gaussian random 
waves have been important in various branches of physics. In particular, a somewhat vague 
heuristic principle due to M.V. Berry [B] asserts that random waves should be a good model 
for quantum chaotic eigenf unctions ( [Zlj contains rigorous results in this direction). Random 
waves in [B] and in much of the physics literature are random Euclidean plane waves of fixed 
energy. For a mathematician, they are defined by putting a Gaussian measure on the infinite 
dimensional space £\ of Euclidean eigenfunctions of fixed eigenvalue A 2 on M m . The Gaussian 
measure satisfies E\ \f\ | 2 = 1, where 1 1 • 1 1 2 is the inner product invariant under the eigenspace 
representation of the Euclidean motion group. 



6 



STEVE ZELDITCH 



On a compact Riemannian manifold, the closet analogue to random Euclidean plane waves 
of fixed frequency is that of random spherical harmonics of fixed degree on S m , where one 
now puts the SO (m + l)-invariant normalized Gaussian measure. In both cases, the measure 
is defined on an eigenspace. For a generic Riemannian manifold, the eigenspaces are of di- 
mension one, so one cannot define an interesting Gaussian measure on the eigenspaces. The 
alternative proposed here and in [Zlj is to replace eigenspaces by the spaces Hn spanned 
by eigenfunctions with asymptotically constant frequency. From the viewpoint of microlo- 
cal (semi-classical) analysis, the analogy is obvious since the spectral projections kernels 
have, to leading order, the same asymptotics as those for spherical harmonics. In |Zlj the 
emphasis was on random orthonormal bases as models for an orthonormal basis of ergodic 
eigenfunctions; while here we only study individual random waves. Riemannian random 
waves of asymptotically fixed energy seem to be a natural global model for random waves on 
a Riemannian manifold without boundary, and the set-up extends naturally to Riemannian 
manifolds with boundary and with fixed boundary conditions on A. When studying local 
behavior, Riemannian waves resemble the Euclidean plane waves of the same frequency. 
More precisely, the scaling limit of Riemannian random waves on length scales A -1 should 
give back the Euclidean plane wave model with eigenvalue 1. This would be the Riemannian 
analogue of the universality result of [BSZlJ. Thus, the natural role of Euclidean plane waves 
seems to be to capture the random behavior of small length scales of order of A~ x . 

Further motivation to study nodal lines of Riemannian random waves has arisen in recent 
conjectures that nodal lines of random two-dimensional Euclidean plane waves of fixed fre- 
quency (and chaotic eigenfunctions) tend to SL e curves [BSt IFGSj IBGS[ 155]. This behavior 
should be sufficiently universal that it should hold for Riemannian random waves of asymp- 
totically fixed frequency on general surfaces. A related conjecture asserts that random nodal 
lines of partial sums of the Gaussian free field (i.e. random Riemannian waves in the cutoff 
ensemble) on a two-dimensional Riemannian surface tend to SLE4 curves [SSJ. Thus, one 
expects different behavior of random Riemannian waves for (almost) fixed frequency and for 
long frequency intervals; although the SLE connection is far outside the scope of this article, 
it does motivate us to consider both ensembles. For recent and deep results on nodal lines of 
random spherical harmonics which are related to conjectures in \BS\ IBGS] we refer to [NSJ . 

The rationale for studying complex zeros of Riemannian random waves is that only in 
the complex domain can we rigorously compare the nodal sets of ergodic eigenfunctions 
and those of random waves. In addition, the complex zeros of complexified Riemannian 
random waves is a higher dimensional generalization of the classical ensembles of of Kac- 
Hammersley of complexified random real polynomials. They may be viewed as sums of 
complexified eigenfunctions on a circle, although the spectral intervals are [0, N] rather than 
[N — 1, N]. Another recent study of complex zeros of complexified real polynomials is the 
thesis of B. MacDonald [Mcj ; however, a significant difference is that the polynomials there 
are orthonormalized in the complex domain rather than the real domain. 

1. Background on densities and correlations of zeros of real Gaussian 

random waves 

In this section, we apply the formalism in |BSZ2j to give explicit formulae for the densities 
of zeros of Riemannian random waves. The same formalism also could be used to give 
formulae for correlations between zeros. 
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1.1. Assumptions on (M,g). We will assume the geodesic flow G l of (M,g) is of one of 
the following two types: 

(1) aperiodic: The Liouville measure of the closed orbits of G t , i.e. the set of vectors 
lying on closed geodesies, is zero; or 

(2) periodic = Zoll: G T = id for some T > 0; henceforth T denotes the minimal period. 
The common Morse index of the T-periodic geodesies will be denoted by (3. 

In the real analytic case, (M, g) is automatically one of these two types, since a positive 
measure of closed geodesies implies that all geodesies are closed. We only need to assume 
(M, g) is real analytic when considering complex zeros. In the C°° case, it is simple to 
construct examples with a positive but not full measure of closed geodesies (e.g. a pimpled 
sphere) . 

The two-term Weyl laws counting eigenvalues of vA are very different in these two cases. 

(1) In the aperiodic case, Ivrii's two term Weyl law states 

N(X) = #{j : \j <\} = c m Vol{M,g) \ m + o(A m_1 ) 

where m = dimM and where c m is a universal constant. 

(2) In the periodic case, the spectrum of vA is a union of eigenvalue clusters Cn of the 
form 

C N = {(^)(iV + J) + fi m , i = l...d N } 

with fi Ni = 0(A^ _1 ). The number d]y of eigenvalues in Cn is a polynomial of degree 
m — 1. 

We refer to [Ho[ ISVj IZlj for background and further discussion. 

1.2. Definition of Riemannian random waves. To define Riemannian random waves, we 
partition the spectrum of a/A^ into certain intervals In of width one and denote by Hj N the 
spectral projections for \/A^ corresponding to the interval In- The choice of the intervals In 
is rather arbitrary for aperiodic (M, g) and as mentioned above we assume In = [N, N + 1]. 
But the choice has to be made carefully for Zoll manifolds. 

In the Zoll case, we center the intervals around the center points + | of the iVth 
cluster Cn- We call call such a choice of intervals a cluster decomposition. We denote by 
djq the number of eigenvalues in I N and put TCn = ranllj^ (the range of Hi N ). Thus, Hn 
consists of linear combinations Yl,j-\ ^i N c j i PN,j of the eigenfunctions {fNj} of a/A^ with 
eigenvalues in In- 

The formalism is simpler in the cutoff ensemble and only requires small modifications from 
the asymptotically fixed frequency ensembles, so we only explain at the end how to modify 
the results in that case. 

1.3. Density of real zeros. The formula for the density of zeros of random elements of 
1~Ln can be derived from the general formalism of |BSZll IBSZ2t IBSZ3] . 

As above, we let \Zf\ denote the Riemannian (m — 1)- volume on Zf. By the general 
formula of [BSZHIHSZ2], 

E\Z fK \=K?{z)dV gi K»(x) = jD(0,£,x)\\£\\d£. (11) 
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We digress to connect this formula with the discussions in |BSZ11 BSZ2, Ne]. In these articles 
||£|| is written a/ det (££*). However, det(££*) = ||£|| 2 in the codimension one case. Indeed, let 
df* be the adjoint map with respect to the inner product g on T X M. Let df x odf* : T X M — ► K. 
be the composition. By det df x o df* is meant the determinant with respect to the inner 
product on T X M\ it clearly equals \df\ 2 in the codimension one case. 
The formulae of [BSZ1[ IBSZ2j (the 'Kac-Rice' formulae) give that 

D(p,bz) = Z n (z)D A (&z), (12) 

where 

Da ^z) = ~ /~r~~~T ex P (13) 
7r m vdet A 

is the Gaussian density with covariance matrix 

A = C - B*A- l B = (C q g , - B q A- l B ql ) , (q = 1, . . . , m) (14) 

and 



In the case at hand, 



Z(.) = = * . (15) 

7rvdetA 7tv A 



A N B N 



(A N ) = E(l 2 ) = V(x,x), 
a N 

(B») t = E(xg,) = ^n,.(.,»)U, 

= E(=E,) = ±^-n & (x, !/ )U, 
g,</ = 1, . . . ,m. 

Making a simple change of variables in the integral (II ip . we have 

Proposition 1.1. [BSZ1J On a real Riemannian manifold of dimension m, the density of 
zeros of a random Riemannian wave is 

K ?M = ,Wn^> J ' J|AWW '^ l|eXPH ^ ))<g - < 16 > 

where A N (x) is a symmetric form on T X M. For the asymptotically fixed freqency ensembles, 
it is given by 

aN ( x ) = T~ \ d x ® d y Ii lN (x,y)\ x=y - — — r -d x U lN (x,y)\ x=y ® d y U lN (x,y)\ x=y ] . 

In the cutoff ensemble the formula is the same except that Hj N is replaced by H.[q,n] ■ 
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2. Zeros of random real Riemannian waves: Proof of Theorem □ 

We begin the proof with the simplest case of the round metric on S m . Throughout this 
article, C m denotes a constant depending only on the dimension. It may change from line 
to line. 

2.1. Random spherical harmonics. To prove Theorem @ on a round S m , we first need 
to evaluate the matrix above when Hn(z,w) is the orthogonal projection onto spherical 
harmonics of degree N. 

Proposition 2.1. Let U N : L 2 (S m ) — > Hn be the orthogonal projection. Then: 

• (A) Tl N (x,x) = Vol \ Sm) d N ; 

• (B) d x U N (x,y)\ x=y = d y Il N (x,y)\ x=y = 0; 

• (C) d x ® d y U N (x,y)\ x=y = mVo ] {sm) X 2 N d N g x . 

Proof. Statement (B) follows from statement (A) since 

= d x IL N (x,x) = d x Il N (x,y)\ x=y + d y Il N (x,y)\ x=y 

and because d x U.N(x, y)\ x=y = d y YlN(x,y)\ x=y . Statement (C) holds because d x <g>d y U.N(x, y)\ x=y = 
CnQx by SO{m + 1) symmetry. To evaluate we use that 

= AIl N (x,x) = 2Trd x ® d y Il N (x,y)\ x=y - 2X 2 N Il N (x,x). 

Here, Trd x ® d y Yl^{x, y)\ x = y denotes the contraction. It equals Ylj=i Wd^NjirfW 2 where 
{^>Nj} is an orthonormal basis. Thus, ttiCn = X 2 n IIn(x,x) and the formula of (C) follows 
from (A). 

□ 

The expected density of random nodal hypersurfaces is given as follows 
Proposition 2.2. In the case of S m , 

K?(x) = C m \ N ~C m N, (17) 
where C m = J Rm |f | exp (-(£, f )) df . 
Proof. By Propositiosn 11.11 and 12.11 we have 

K?{x) = f \\A N (x)^\\ex P (-(t,0)dti, (18) 

where 



□ 
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2.2. Random Riemannian waves: proof of Theorem [H We now generalize the result 
to any compact C°° Riemannian manifold (M, g) which is either aperiodic or Zoll. As in the 
case of S m , the key issue is the asymptotic behavior of derivatives of the spectral projections 

U lN (x,y)= ^{x)^M- ( 19 ) 

j:\jeI N 

Proposition 2.3. Assume (M,g) is either aperiodic and I N = [N, N + 1] or Zoll and I N 
is a cluster decomposition. Let Hj N : L 2 (M) — > H.^ be the orthogonal projection. Then: 

• (A) Ui N (x,x) = Vol{ M tg)) d N (l + o(l)); 

• (B) d x U lN (x,y)\ x=y = d y Il N (x,y)\ x=y = o(N m ); 

• (C) d x ®d y U lN (x,y)\ x=y = Vol ^ g)) \%d N g x {l + o{l)). 

In the aperiodic case, 

(1) U m (x,x) = C m A m + o(A m - 1 ); 

(2) d x ® d y U [0 , X ]{x,y)\ x=y = C m X m+2 g x + o(A m+1 ). 

In the Zoll case, one adds the complete asymptotic expansions for Yii N over the N clusters 
to obtain expansions for II ^ . 

Proof. Asymptotic formulae of type (A) are standard in spectral asymptotics, and we refer 
to |DG[ IHo] for background. Asymptotics of type (C) were worked out in [Z2j (Theorem 
2) in the special case of a Zoll metric. However, much of the calculation goes through for 
any compact Riemannian manifold. It does not appear however that (B) has been stated 
before or that (C) has been previously discussed on general Riemannian manifolds, although 
the techniques are standard. These asymptotics are dual to the heat kernel asymptotics in 
[BBGJ, but they are sharper because we are using spectral intervals for \/A of fixed width 
rather than intervals of the form [0, A], which are dual to heat kernel asymptotics. Thus, 
we need two term asymptotics for long intervals in order to obtain asymptotics on short 
intervals. 

We follow the standard Tauberian method of |DGj . Proposition 2.1, for studying the 
spectral asymptotics. We consider the spectral measure 

d x ii[o,\](x,y) = x^ipxji^fXjiy), (20) 

whose integral over the interval In equals Ui N (x,y), and the derived measures 
(a) dxU[ ,x](x,x) = ^2j5(X - \j)yx 3 {x) 2 

< (b) d x d x Il[ 0; x]{x,x) = 2^2 j 8(X- \j)ip Xj {x)d(px j {x) . (21) 

(c) dxd x <^d y U[o ) x](x,y)\ x = y = J2j^-^j)dipx j {x)<^id(px j (x). 

We now introduce a cutoff function p € «S(R) with p 6 supported in sufficiently small 
neighborhood of 0. We also assume p = 1 in a smaller neighborhood of 0. Then there exists 
an expansion in inverse powers of A with coefficients smooth in (x,y): 

00 

p* dxU [0 ,x](x,x) = ^p(A - Xj)(p\.{x) ~Y,W*- x - k u k {x), (22) 

3 k=0 

where u>k are smooth in x, and u = 1- 
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We briefly recall the proof of ( 1221) and of (A): We have 

p*d x U [0>x] {x,y) = [ p(t)e itx U(t,x,y)dt, (23) 

where U(t, x, y) is the Schwartz kernel of the wave group U(t) = e~ itVE . We use a small-time 
parametrix for U (t, x, y) near the diagonal of the form 

U(t,x,y)= I e~ itl ^y e i{ii ' exp y 1{x)) A(t,x,y,£)d£ (24) 

J T*M 

where \£\ 9x is the metric norm function at x, and where A(t, x, y, £) is a polyhomogeneous 
amplitude of order which is supported near the diagonal. Setting x = y gives 

p*d x U [0tX] {x,x) = / p(t)e itx e- im syA{t,x,x,£)d£dt. (25) 
Jr Jt*m 

As in |DGj . we pass to polar coordinates r = \£\ g , change variables 9 —>■ X9 and apply the 
stationary phase method to the drdt integral to obtain (j22"j) . 

For (C), we apply the method of [Z2] (see (3.6) - (3.7)). We denote the phase of U(t, x, y) 

by 

<fi(t, x ,y,£) = (^expy 1 ^)) -t\€\gy 
Then applying d x <S> d y \ x=y to the integral produces a highest order term given by a universal 
constant times 

d x ip(t, x, y, £) ® rf y p(t, x, y, 1^ = £ ® ^, 

since a (t,a;, x,£) = 1 (cf. [DG]). It also produces lower order terms (i.e. of order < 1) in 
which at least one derivative falls on the amplitude. If we the put the d£- integral in polar 
coordinates £ = ru, we obtain an expansion 

oo 

p*d x d x ® d y U m (x, y)\ x=y ~ A m+1 " fc 5 fc (x), (26) 

k=0 

with the leading coefficient 

B o( x ) = T/ uZ \ [ u®udp x (uj) = jtj- g x . (27) 
Vol{M,g) J 3 . M Vol{M,g) 

Here, dp x is the Euclidean area element induced by g on S*M. 

We now draw the conclusions for (A) - (C). In the Zoll case, (A) and (C) are already 
proved in detail in [Z2J. In the Zoll case, there exist complete asympotic expansions for the 
spectral sums and we may deduce (B) as well from the smoothed expansion by a modification 
of the proof of |Z2j . Theorem 2. In the aperiodic case, the measures (A) and (C) are positive 
and we may apply the Fourier Tauberian theorems of |Hol ISVj (see the Appendix in $6]), or 
alternatively the remainder estimate of Ivrii, to obtain two-term expansions: 

(1) U m (x,x) = C m \ m + o{\ m - 1 )- 

(2) d x <g> d y U [0 , x] (x,y)\ x=y = C m X m+2 g x + o(X m+1 ). 

We then subtract the expansions across the interval 7jy to obtain the stated result. We 
note that the drop in degree is encoded in dx ~ iV m_1 . 
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We still need to analyze (B). Since (B) equals d x Ili N (x, x), we could take the derivative 
in x of (I25I) to obtain 

p*d x d x IL [0>x] (x,x) = ( [ p(t)e itx e- u ^yd x A(t,x,x^)d£dt (28) 

Jr Jt*m 

We then have an expansion similar to that of (j22p except that the amplitude is now the 
one-form d x A. However, the Tauberian theorems do not apply to (B) since d\d x Il[o^x\(x, x) 
is not a positive measure. In the Zoll case, we have a complete asymptotic expansion of (A) 
and its x-derivative gives that of (B). But on a general Riemannian manifold one cannot use 
this approach. 

Henceforth we assume (M,g) is aperiodic. In this case, we use the fact that, for each k, 
B k = t|^-A -1//2 is a bounded pseudo-differential operator and 

d 

—U [0tX] (x,y)\ x=y = 2^ (B k <p x .(x)) <p\ s (x). 

j:\j<\ 

We write 

2(B k <p x .(x))ip Xj (x) = ((I + B k )<p x .(x)) 2 - (B k <p Xj (x)) 2 - <p 2 Xj . 

We then substitute the right side into the summatory function in Xj to obtain three 
asymptotic expansions to which the Tauberian theorems apply. We start with 



p * d x ((I + B k ) x <g> (I + B k ) y V AII[o,a] (x, 

= ((/ + B k ) x + B k ) y J M p(t)e- itx VAU (t, x, y) \ x=y dt 



((/ + B k ) x ® (I + B k ) y J R p (t)e itx £ t U(t, x, y) \ 



x=y 



{{I + B k ) x ® (/ + B k ) y ( J R f T * M (ip'(t) - Xfie^e-Wnj&^WAit, x, x, £)d£dt 



I x=y 

(29) 

To leading order, application of ((J + B k ) x ® (/ + Sfc)j/ under the integration sign multiplies 
the leading term of the amplitude by (1 + b k )(x,d(p t ) where <p t is the phase (cf. e.g. the 
'fundamental asymptotic expansion' of [T]). We then apply the stationary phase method 
and due to the extra factor of A in the amplitude obtain 



2:,:a • a A, ((/ + B k )^p Xj {x)f = C m A m+1 f s$M (l + b k (x,u)fdp(uj) + o(A m ). (30) 

The o(A m ) remainder holds as in the scalar case because the geodesic flow is aperiodic 
[DGl IHol ISVj . We then repeat the calculation for B k and for / and subtract. We clearly 
cancel the leading term, leaving the remainder o(A m ). When we subtract the interval [0, N] 
from [0, iV + 1] we obtain o{N m ). 

□ 



2.3. Proof of Theorem The generalization of Proposition 12.21 to a general Riemannian 
manifold is the following: 
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Proposition 2.4. For the asymptotically fixed frequency ensemble, and for any C°° (M,g) 
which is either Zoll or aperiodic ( and with In as in Proposition \2.3\) , we have 

K?{x) = ^x^/ K Jlel|exp(-^(e,0)^ + o(l) 

V J (31) 

~ C m N, 

where C m = J Rm ||£|| exp (— (£, £)) d£. The same formula holds for the cutoff ensemble. 

Proof. Both on a sphere S m or on a more general (M, g) which is either Zoll or aperiodic, 
we have by Propositions 12.11 resp. 12.31 and the general formula for in §1.31 that 



A N ( s 1 ( (1 + 0(1)) 0(1) 

[Z) ~ Vol(M,g)\ o(l) N 2 ^(l + o(l)) 



(32) 



(34) 



It follows that 

a n = c n_ B N* {A Nyi B N = vd * M N 2 g x + o(N). (33) 

Thus, we have 

K?W ~ ^W^L m \\A N (x) 1/2 t\\ew(-(t,0)dZ 

= £ J Rm + o(l))(x)m I exp (-<£, 0) de, 

where o(l) denotes a matrix whose norm is o(l). The integral tends to L m 1 1£| | exp (— (£, £)) 
as iV — » oo, completing the proof. 

□ 

So far, we have only determined the expected values of the nodal hypersurface measures. 
To complete the proof of Theorem [Tj, we need to prove: 

Proposition 2.5. If(M,g) is real analytic, then the variance of j^X^ is bounded. 

Proof. By a theorem of Donnelly- Fefferman [DF] , for real analytic (M, g) , 

Cl A < U m -\Z VX ) < C7 2 A, (A^ A = AVa; c 1s C7 2 > 0). (35) 

Hence for any f^ G Tir^, —Zf N has bounded mass. Hence, the random variable j-X^ is 
bounded, and therefore so is its variance. 

□ 

Remark: 

The variance of -^X^ is given by 

Var(±X») = -L [ [ (K»{x,y) - K?{x)K»{y)) ^(x)^(y)dV g (x)dV g (y), (36) 
*n a n j m j m 
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where K^{x,y) = E 7jv (Z/ JV (x) (g> Zf N (y)) is the pair correlation function for zeros. Hence, 
boundedness would follow from 

-V/ / K»(x,y)dV g (x)dV g (y)<C. (37) 

There is a formula similar to that for the density in Proposition 11.11 for K^ix^y) and it is 
likely that it could be used to prove boundedness of the variance for any C°° Riemannian 
manifold. But we leave this for the future. In the Kahler case, asymptotic formulae for the 
variance of smooth linear statistics are given in [SZ2, ISZ3J . but the method does not apply 
in the real case. 



2.4. Random sequences and proof of Corollary [21 We recall that the set of ran- 
dom sequences of Riemannian waves of increasing frequency is the probability space Ti^ = 
UzvLi Win w hh the measure = []^ =1 7^. An element in TL^ will be denoted f = {Jn}- 
We have, 

\(^z fN ,n<Y 

By a density argument it suffices to prove that the linear statistics j^(Zf N , if)) — y^pj^y f M ipdV g 
almost surely in Tioo. From Theorem [TJ we have: 

Corollary 2.6. (%) lim^oc £ E k <N E (ir k x i) = vsm JmW; 
(it) Var(j^Xff) is bounded on Ti^. 

Since jj^X^ for {,N = 1,2, . . .} is a sequence of independent random variables in 7Yoo 
with bounded variances, the Kolmogorov strong law of large numbers implies that 



almost surely. 



3. Analytic continuation of eigenfunctions and spectral projections 

We now turn to complex zeros of analytic continuations of random Riemannian waves. 
Before getting into the statistics of zeros, we need to recall the basic results on analytic 
continuation of eigenfunctions and to introduce the basic two-point kernels. 

As mentioned in the introduction, for each analytic metric g there exists a unique plurisub- 
harmonic exhaustion function p on Mq inducing a Kahler metric u g = iddp which agrees 
with g along M. We recall y/p(C) = ^ r c(C C) where r(x, y) is the distance function and r c is 
its holomorphic extension to a small neighborhood of the anti-diagonal ((, () in M<c x Mq, Jp~ 
is a solution of the homogeneous complex Monge- Ampere equation {dd^fp) m = away from 
the real points. We refer to |GS1[ IGS2[ ILSl} IGLS} |Z3] for further background on Grauert 
tubes and adapted complex structures on cotangent bundles of analytic Riemannian mani- 
folds. 

The eigenfunctions (pj admit holomorphic extensions tpj to the maximal Grauert tube 
|Bou|IGS2] . As a result, one can holomorphically extend the spectral measures dLTro^]^, y) = 
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^■5(A — Xj)(fij(x)(fij(y) of V~K. The complexified diagonal spectral projections measure is 
defined by 

dxnf 0iX] (c,() = J2 6 ( x - x M(O\ 2 - (38) 

j 

Henceforth, we generally omit the superscript and write the kernel as (). This kernel 

is not a tempered distribution due to the exponential growth of |^(C)| 2 - Since many as- 
ymptotic techniques assume spectral functions are of polynomial growth, we simultaneously 
consider the damped spectral projections measure 

dxP{oA)(CO = ^K^IftfCOl 2 , (39) 

j 

which is a temperate distribution as long as y / p(C) < r - When we set r = y / p(C) we omit 
the r and put 

cUWCC) = - ^)e" 2 ^ (C)Aj |^ C (C)r. (40) 

i 

The integral of the spectral measure over an interval / gives ILi(x, y) = X^-a ei Vji^Pjiv)- 
Its complexification gives the kernel ([9]) along the diagonal, 

n,(C,0 = £ I^OI 2 , (4i) 

and the integral of fl39|) gives its temperate version 

^ r (c,o= e e ~ 2rAj M(or, (42) 

or in the crucial case of r = v /p(C), 

^/(C,C)= £ e- 2 ^|^(C)| 2 , (43) 

3.1. Poisson operator as a complex Fourier integral operator. The damped spectral 
projection measure d\ P^ x ^((,() ( |39i) is dual under the real Fourier transform in the t 
variable to the restriction 

U(t + 2ir,C0 = Y^^ 2T+U)X3 H(0\ 2 (44) 
j 

to the anti-diagonal of the mixed Poisson-wave group. We recall that the kernel U (t+ir, x, y) 
of the the Poisson-wave operator e *(*+ lT )v / A ac i m its a holomorphic extension in the x variable 
to the closed Grauert tube M T . The adjoint of the Poisson kernel U(ir,x,y) also admits an 
anti-holomorphic extension in the y variable. The sum (jUJ) are the diagonal values of the 
complexified wave kernel 

U(t + 2iT,(,(') = J M U(t + iT,(,y)E(iT,y,(')dV g (x) 

(45) 

= E i e<- 9r+ *^^(C)^(C')- 
We obtain fj4"5]) by orthogonality of the real eigenfunctions on M. 
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Since U (t + 2ir, (, y) takes its values in the CR holomorphic functions on dM T , we consider 
the Sobolev spaces O sA ~{dM T ) of CR holomorphic functions on the boundaries of the 
strictly pseudo-convex domaisn M £ , i.e. 

o s+ ^{dM T ) = w s+I ^(dM T ) n o(dM T ), 

where W s is the sth Sobolev space and where 0(dM e ) is the space of boundary values 
of holomorphic functions. The inner product on O (dM T ) is with respect to the Liouville 
measure 

dfir = {idd^p) m ~ l A d c Jp. (46) 
We then regard U(t + ir,(,y) as the kernel of an operator from L 2 (M) — > O (dM T ). It 
equals its composition Il r o U(t + ir) with the Szego projector 

n r : L 2 {8M T ) -> O {dM T ) 

for the tube M T , i.e. the orthogonal projection onto boundary values of holomorphic func- 
tions in the tube. 

This is a useful expression for the complexified wave kernel, because Il r is a complex 
Fourier integral operator with a small wave front relation. More precisely, the real points 
of its canonical relation form the graph As of the identity map on the symplectic one 
S r C T*dM T spanned by the real one-form d c p, i.e. 

E r = {(C; rd°p(()), C e dM T , r > 0} C T*(dM T ). (47) 

We note that for each r, there exists a symplectic equivalence S r ~ T*M by the map 
(£, rd c p(Q) — > (E^ 1 ^), rat), where a = £ ■ dx is the action form (cf. |GS2]). 

The following result was first stated by Boutet de Monvel (for more details, see also 
[GS21IZ5]). 

Theorem 3.1. |Boul IGS2j n e o U(ie) : L 2 (M) — * 0(dM £ ) is a complex Fourier integral 
operator of order associated to the canonical relation 

T = {(y,r ] ,L £ (y,r ] )}cT*Mxi: £ . 

Moreover, for any s, 

Tl £ o U(ie) : W S (M) -> O s+m ^{dM e ) 
is a continuous isomorphism. 

We obtain the holomorphic extension of the eigenfunctions ip\ of eigenvalue A 2 by applying 
the complex Fourier integral operator U(ir): 

U{ir)ip x = e-'VS- (48) 
Combining Theorem 13. II with a stantard Sobolev estimate, we obtain: 

Corollary 3.2. [Boy] IGLSj Each eigenfunction (fx has a holomorphic extension to B*M 
satisfying 

sup |^ A (C)|<C £ A m+ V A . 

The order of magnitude reflects the fact that E{ir) smooths to order — since it is a 
bounded operator from a manifold M of real dimension m to one dM T of dimension 2m — 1. 
We will also need the following Lemma from [Z3j (Lemma 3.1): 



REAL AND COMPLEX ZEROS OF RIEMANNIAN RANDOM WAVES 17 

Lemma 3.3. Let a e S°(T*M - 0). Then for all < e < e , we have: 

U(te)*U £ aU £ U(ie) G ty~ m * 1 (M), 

with principal symbol equal to a(x,£) \£\ g 2 . 

Here, \1/ S (M) is the class of pseudo-differential operators of order s and Eq is the radius of 
the maximal Grauert tube. 

4. Complex zeros of random waves: Proof of Theorem 

Our first result is a formula for the expected distribution of complex zeros. We retain the 
same ensembles and notation from the real case. 



Proposition 4.1. For any test form ip , 

j(E JN [z f c},i;) = i<09io g ii/ w (c,c),v> 

Proof. By the Poincare-Lelong formula Zjc = ^ddlog |/ C (C)| 2 ; hence we have 

V-yAZf] = ^f lN dd\og |/ C (C)| 2 = ^ddE, N log |/ C (C)| 2 . (49) 

We write / = Y.\ j( n N c Nj^Nj and then complexify to obtain / C (C) = Yjx^In c Nj<PNj(0- 
Thus, we need to calculate E 7iV log | ^2 Xj ei N c NjfNj\- 
Define 

$ N = (<p N1 ,...,<p Ndrr ) :M^C dN 
to be the vector of eigenf unctions from the orthonormal basis of eigenf unctions in the spectral 
interval. Then, 

d N 

l^(C)l 2 = £l^i(C)l 2 = n? w (C,C)- 

We use the notation: 

= [/(c) + iv((), ; u(0, v(0 e R d », \u\ 2 + \v\* = i. (so) 

Lemma 4.2. E 7iV log |/ C (C)| 2 = log 11/^^,0 + G N ((,(), where G N ((,() is the uniformly 
bounded sequence of continuous functions on Mc given by 

G N ((, C) = r'(i) + r(|) logmax{||?7 + JV\\\ \ \U - JV|| 2 }, 
where J is the natural complex structure on the real 2-plane spanned by U, V in M. dN . 

Proof. We have, 

E 7iv log | Exjei N cnjVnj I = f R * N lo S I ( c ' ®n) \dv k ( a ) 

= log|<^(C)| 2 + GV(C,C), 

where 

Gn((,() = We-l^logKcHg)!^ (51) 
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Hence it suffices to calculate G n ((,()■ 

The span of U, V is a real two-dimensional plane in M. dN when U, V are linearly independent; 
we will consider the other case in the remark below. We rotate coordinates in the integral 
for each £ so that 

U — (cos 9 coscp) ei, V = (cos0sin<£>)ei + sin#e 2 , 

where {ej} is the standard basis of lR m and where 9, tp are real angles depending on £ and 
fc. Then, 

(a, i — ;ftttt) = ctifcos coscp) + i cos 6 1 sin 09) + ia 2 sin 0. 
K(C)I 

After rotating coordinates in the integral, the integrand depends only on 0,1,0,2, so we may 
integrate out the remaining variables 03, ... , Od N and obtain 

Gn(CX) = / M2 e-H 2 /2l g|( a; ||^)| da 

*£r(C) 



io°° e 2 Jsi l °g\( ru ']^o\)\ rdrduJ 

r'(|) + r(|) f 51 log |§§)|^ 



r'(f) + r(±)<^(c,c), 



>2/ 1 \2> 

where 

G2 ((, — Jci log |^i (cos 6 1 cos </? + i cos 6* sin </?) + iu 2 sin 0|da; 

= J 27r log \ae iu + (3e~ iu \du = log \ae i2u + (3\du (52) 
= max {log |a|, log|/3|} 

with 

2a = (U + iV, ei - ze 2 ) = cos #e i</3 + sin 9, 2(3 = (U + iV, ei + ie 2 ) = cos 9e i{p - sin 5. 

In the two dimensional real space spanned by e\, e 2 we may define a complex structure by 
Je\ = e 2 , Je 2 = —e\ and then 

G^(CC) = logmax{||f/ + JK|| 2 , ||C7 - jy|| 2 }. 



For any two vectors U, V with ||t/|| + 

1 < ma,x{\\U + JV\\ 2 , \\U- JV\\ 2 } = max{l + 2(U, JV), 1 - 2{U, JV)} < 3. 

Therefore the logarithm is bounded above and below and is clearly continuous. This com- 
pletes the proof of Lemma 14.21 □ 

To complete the proof of Proposition 14.11 it suffices to integrate the dd by parts onto / in 
the Gn term, and use that Gat is uniformly bounded. □ 

Remark: When (eMwe have V = and \\U\\ = 1 so the calculation is slightly different. 
It is possible that there exist other ( 6 Mc such that U, V are linearly dependent. When 
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<j>c 

V = 0, we put U = C\ and obtain the limiting case when <p = 0, (a, j^cf^) — a i- We integrate 
out all but ai and get 

Gn(CX) = J M e-H 2 /2l g| G | rfa , CgM . 

For sufficiently Grauert tubes, it is impossible that U — 0, but our analysis has not ruled 
out that U, V may be linearly dependent, i.e. the case 6 = 0. But the calculation above goes 
through without change in this case. 

5. Complexified spectral projections and wave group 
In view of Proposition 14. 1[ the remaining step in the proof of Theorem [3] is to show: 
Proposition 5.1. 

ilo g n [0lA] (C,C)- Vp, (53) 

and 

ilo g n Jjv (C,C)^ sft. (54) 

For both the long [0, A] and short In intervals, the upper bound is rather simple to prove. 
The lower bound is more difficult, primarily for the short intervals. As discussed above, in the 
Zoll case, the eigenvalues lie in clusters of width N~ l around an arithmetic progression, so 
unless the intervals In are centered along this progression there need not exist any eigenvalues 
in In and the logarithm would equal — oo. We have by definition centered the intervals on 
the arithmetic progression, but clearly this choice must play a role in the proof. 

5.1. Proof of Proposition l5TTl We first give a simple comparison between the logarithmic 
asymptotics of Ui N to those of Pj N (cf. ( 1431) . 

Lemma 5.2. For any t = y/p(() > 0, there exists C, c > so that 

j + 2 y/p(c) + ~ log p [XM1] (c, C) < \ log n [A)A+1] (c, C) < 2v/p(C) + \ log P[a,a +1] (C, C) + p 

hence 

piogn Jjv (c,C) = 2Vp(C) + pio g P 7jv (C,C) + o(±), 

where the remainder is uniform. 

Proof. Since — 1 < N — \ N j < 1 for \ N j G In, w e have 

d N d N 

^v^ge- 2 ^,!^ (c)|2 < UiAc - c) < Ce^^Y,e-^^%{Q\\ (55) 



Next we observe that the upper bound in both the cutoff ensemble 



□ 



lim sup i log n [0iA ] (C, C) < \J~P (56) 

follows immediately from Corollary 13.21 and the fact (cf. Lemma 1572]) that II^a] < e X ^P{o,x\- 
It is similar but simpler in the fixed frequency ensemble. 
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Hence, it remains to prove is the lower bounds. As mentioned above, the global nature 
of the geodesic flow must play a role in the proof. The route we choose uses the minimal 
amount of information necessary to obtain the result: we only use the behavior of the spectral 
projections in the real domain to deduce a lower bound in the complex domain. We follow 
the method of [SZlj . 

Put V/v(C) := 77 l°gn[jv,jv+i](C) 0- Then Vn is plurisubharmonic in M £ , and we would 
like to show that Vn — > ^fp in L 1 (M e ). If not, we can find a subsequence {V/v fe } with 
ll^Vfc — \fP\^{M e ) > 5 > 0. Since {Vn} is uniformly bounded above, a standard result 
on subharmonic functions (see [Ho[ Theorem 4.1.9]) implies that sequence {V^} either 
converges uniformly to — oo on M £ or else has a subsequence which is converg ent in L>{M £ ). 
The first possibility cannot occur since Vjv(x) — > on the totally real submanifold M. 

Therefore there exists a subsequence, which we continue to denote by {Vat*.}, which con- 
verges in L 1 (M e ) to some G G L l (M £ ). By passing to a further subsequence, we may assume 
that Vjy fc — > G converges pointwise a.e. in M £ . Let V(Q = lim sup^^ V^ k (z) (a.e) and 
let 

V*(z) : = lim sup V(w) 

w—*z 

be the upper-semicontinuous regularization of V. Then V* is plurisubharmonic on M e , 
V* < yfp (a.e.) and V* = V a.e. 

Since \\VN k — ^/pWl 1 ^ 1 ) > 5 > 0, we know that V* ^ J~p. Hence, for some e > 0, the open 
set U £ = {z G M £ : V* < ^ — s} is non-empty. Let U" be a non-empty, relatively compact, 
open subset of U £ . Then by Hartogs' Lemma, there exists a positive integer K such that for 
C G U", k>K,V Nk <^/p- e/2, i.e. 

n^nK.C^e^-*, (eu", k>K. (57) 

By Lemma 15.21 it follows that 

P N (C0<e- ENk , Cell", k>K. (58) 

We now let x C^°(U") be a smooth bump function which equals one on some ball 
B C U". By (EHD, we have 

/ X P N ((XW<e-^. (59) 

JdM £ 

We observe that 

= E j: A je [iV,7V+l](X^(^V / P(C)) ( ^A j) ^(^v / P(C))^A J )L2(OM £ ) 

(60) 

= E j: A je [V,iV+l]( f/ ( i V / P(0)*X^(«V / P(0)^A J ^A J )L2 (M) 

= Trn % ^(^Vp(c))*xf/(^ v / p(C))- 

By Theorem 13.11 and Lemma I3TB1 f^(^y / p(C))*X^(*v / ^(0) i s a pseudo-differential operator 
with principal symbol equal to |£| 2~x(^, £)■ For intervals In satisfying our assumptions, 
it follows from a standard local Weyl law [Ho] (see also [Zlj for further discussion in the 
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present context) that for a zeroth order pseudodifferential operator A, and for a Laplacian 
which is either aperiodic or Zoll, 

TrU lN A = C m N m - 1 [ xdli + o^N™- 1 ). (61) 

J B*M 

n-l 

Putting A = y/A U(t^p(C)YxU^(0) 

TrIL Ijf U(iy/p(Q)* X U{iJp(£)) ~ OJV"^ [ x dfi. (62) 

Jb*m 

The latter asymptotics contradict the exponential decay of (l59l) . 

□ 

6. Appendix on Tauberian Theorems 

We record here the statements of the Tauberian theorems that we use in the article. Our 
main reference is |SVj . Appendix B and we follow their notation. 

We denote by F + the class of real- valued, monotone nondecreasing functions N(X) of poly- 
nomial growth supported on R + . The following Tauberian theorem uses only the singularity 
at t — of dN to obtain a one term asymptotic of N(X) as A — ► oo: 

Theorem 6.1. Let N e F + and let ip e <S(M) satisfy the conditions: ip is even, ip(X) > 
for aliAeK,^ , and ^(0) = 1. Then, 

^*dN{X) < AX U \N{X)~N*tfj{X)\<CAX u , 

where C is independent of A, X. 
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